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We study theoretically the phase diagram of perovskite 
manganites taking into account the double degeneracy of the 
eg orbitals in a Mn^"*' ion. A rich phase diagram is obtained 
in the mean field theory at zero temperature as functions of 
X (hole concentration) and Jg (antiferromagnetic interaction 
between t2g spins). The global features of the phase diagram 
is understood in terms of the superexchange and double ex- 
change interactions, which are strongly depends on types of 
the occupied eg orbitals. The strong electron correlation in- 
duces the orbital polarization, which controls the dimension of 
the conduction band. A sequential change of the spin and or- 
bital structures with doping holes is consistent with the recent 
experiments. In particular, metallic A-type (layered) antifer- 
romagnetic state is found for x ~ 0.5 with the uniform d^2_y2 
orbital ordering. In this phase, calculated results suggest the 
two-dimensional conduction and absence of the spin canting, 
which are observed experimentally. Effects of the Jahn- Teller 
distortion are also studied. 

71.27.-Fa, 75.30.-m, 75.30.Et 



I. INTRODUCTION 

Doped manganites Ri-xAxMnOa (R=La, Pr, Nd, Sm 
; A= Ca, Sr, Ba) have recently attracted considerable in- 
terests from the viewpoint of a close comifiCjiion between 
the magnetism and the electric transport .Eru Theoretical 
studies of the double epcghange interaction have been de- 
veloped long time agoO~El and explained the, emergence of 
the ferromagnetism in doped manganites.u However re- 
cent systematic experimental studies have revealed more 
rich phase diagrams in this system. 

FIG. 1. 

Observed phase diagrams of Lai-ajSr^^MnOs (wider 
band width system) and of Pri-ajSr^jMnOa and 
Ndi_2:Sra;MnOjh|-|(|ii&rrower band width systems) are 
shown in Fig|^.HiJ'll2l In the parent compound (x = 0.0), 
the layered antiferromagnetism (spin A-type AF) ac- 
companied|-with a distortion in the MnOe octahedron is 
realized.Oila By moderately doping holes, the insulator 
with spin A-type AF changes into an insulator with a fer- 
romagnetism (spin F-type) around-E = 0.125, and to a 
spin F-type metal at a; - 0.1751^^ In Pri^^Sr^MnOa 
and Ndi-xSr^jMnOa, a metallic phase with snirL_A-type 
AF is found recently to appear for x > 0.48Bt3'lij With 
further increasing of a; (a; ^ 0.6), rod type antifcrromajg^ 
netism (spin C-type AF) is found in Lai-xCajjMnOsEJ 



and Ndi_a;Sr^MnO30 Finally, for x = 1.0, the three di- 
mensional AF (spin G-type AF) is realized. In addition, 
in the narrower-band systems, the charge ordered state 
accompanied by the spin ordering is recognizedjajsar the 
commensurate value of x (x=0.5, 0.75 etc.),t£HlZl where 
the orbital is also supposed to be ordered. 

In order to reveal the origin of the unique magneto- 
transport in this system, it is essential to understand 
the above rich phase diagram. However, it cannot be 
explained by the conventional scenario based on the dou- 
ble exchange interaction. This discrepancy should be at- 
tributed to ingredients neglected in the conventional the- 
ory, namely the anisotropic transfer intensity originated 
from the eg orbital degrees of freedom, electron-electron 
interactions, and the electron-lattice interaction (Jahn- 
Teller (JT) effect). Especially, in the narrower-band 
systems, the kinetic energy of Cg carriers is suppressed 
and the above interactions become more important. At 
X = 0.0, the spin and orbital structures have been stud- 
ied theoretically by taking into account the above inter- 
actions sin ce , t he pioneering works by Goodenough and 
Kanamori£3il3 However, when we focus on the origin of 
the spin structure at x = 0, i.e., the A-type AF, and roles 
of the JT distortion on it, situation is still controversial. 
In one side of the theoretical investigation, the spin or- 
dering is attributed to the strong electron-electeoaJnter- 
action and the doubly degenerate eg orbitalsEMcS The 
ferromagnetic superexchange interaction, which is neces- 
sary to explain the spin alignment in the afo-plane, is orig- 
inated from the degenerate orbitals. and—the Hund cou- 
pling interaction between them.E3rJ'EaE3 Even without 
JT effect, the A-ytpe spin alignment is derived by the su- 
perexchange interaction and the anisotropy of the trans- 
fer integral due to the orbital ordering. However, a type 
of the orbital ordering theoretically obtained disagrees 
with that expected from a type of the JT distortion, i.e., 
dsx^-r^ and d^y2_j.2. The another side of the thcorcti^ 
cal investigation is based on the Hartree-Fock theorjtZl 
and the first-principle baiid^ calculation base on the lo- 
cal density approximation£3 where the JT distortion is 
indispensable to reproduce the observed spin structure 
through introduction of the orbital polarization. With- 
out the JT distortion, the system becomes a ferromag- 
netic metal. 

As for the doped case(a; ^ 0), the situation is even 
more controversial. An issue of main interest here is the 
origin of the colossal magnetoresistance (CMR) observed 
near the ferromagnetic transition temperature Tc. For 
this purpose identification of the dominajit interactions 
is of primary importance. Millis et alE3 attributed an 
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insulating behavior of the resistivity p{T) above to 
formation of the small JT polaron. They assumed the 
strong Hund's coupling, but the other Coulomb interac- 
tions were neglected. The characteristic JT interaction 
is about leV which is comparable or smaller than the 
band width. Then it is reasonable that the small polaron 
formation disappears below Tc and p{T) shows metallic 
behavior. However we argue that this picture does not 
explain the following anomalous features observed in the 
low temperature ferromagnetic state (T « Tc): (i) p{T) 
is fitted by p{T) = po + AT^ with the coefficient A being 
large, i.e., A ^ 500 p^l cm/ K^a (ii) the optical conduc- 
tivity criuj) is domin|ated by the incoherent part with a 
small Drude weight (iii) photo-emission experiments 
show only a small discontinuity at the Fermi edge.EJ We 
regard these features as signatures of the strong correla- 
tion in the doped Mott insulator and the Coulomb inter- 
actions remain strong even in the metallic state. Con- 
sidering the strong Hund's coupling, which causes the 
perfect spin polarization, it is reasonable to assume the 
strong correlation between Cg electrons. In the ferromag- 
netic state, the orbital degrees of freedom play a similar 
role to that of spins in the usual doped Mott insulator. It 
is the additional ingredient important in the manganites. 

In this paper we study a phase diagram of perovskite 
manganese oxides. The double degeneracy of the eg or- 
bital, the anisotropy of the transfer integral based on it 
and the strong electron correlation are taken into account 
in the model. The spin and orbital phase diagram at zero 
temperature is obtained by the mean field approximation. 
The sequential change of the spin structure with doping 
of holes, I.e., A^F^A^C^G, is well reproduced 
by the calculation. We found that the strong Coulomb 
interactions experimentally suggested induce the perfect 
orbital polarization which plays an essential role to de- 
termine the spin structure. The orbital structure is also 
changed with doping of holes and controls the dimension- 
ality of the conduction bands. In particular, the metallic 
spin A- type AF phase is found around a; '--^ 0.5 where 
the orbitals are aligned as dx^-.y2. In this phase, the in- 
terlayer electron transfer is forbidden both by the spin 
and orbital structures and the spin canting is absent. 
Both theoretical ucedictions are consistent with the re- 
cent experiments.S'Ej Roles of the JT distortion in the 
undoped case is also studied. 

In Sec. II, we introduce the model Hamiltonian and 
the formulation of the mean field calculation. Results 
of the numerical calculation are presented in Sec. III. 
Sect. IV is devoted to discussion and conclusion includ- 
ing comparison with the previous works short version 
of this paper has been already published,E3 but this paper 
contains additional and more detailed results. 



II. MODEL AND FORMULATION 
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where Hk is the kinetic energy of Cg electrons, Hnund is 
the Hund coupling between Cg and t2g spins, and Hon- site 
represents the on-site Coulomb interactions between Cg 
electrons. t2g spins are treated as the localized spin with 
S = 3/2. The AF coupling between nearest neighboring 
t2g spins is introduced in Hs to reproduce the AF spin 
ordering observed at x = 1.0. Using an operator dj^^ 
which creates an eg electron with spin a (= t, i) in the 
orbital 7 (= a{d^2 _y2) , b[d3^2 _^2)) at site i, each term of 
Eq.(n]) is given by 
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tjj in Hk is the electron transfer intensity between near- 
est neighboring sites and it depends on kind otpccupied 
orbital and the direction of a bond as follows H 
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(5) 



(6) 



(7) 



to is the electron transfer intensity between d-^^^^j.^ or- 
bitals along the z-direction. The spin operator for the Cg 
electron is defined as Se^i = 5 '^-ya^apd-^p with the 

7a/3 

Pauli matrices (Ja/s- St2gi denotes the localized t2g spin 
on i-site with S = 3/2. The last term in the Hamiltonian 
Hon-site cousists of the following three contributions: 



Ho 



— Hjj + Hw + Hj 



(8) 



where Hu and Hu' are the intra- and the inter-orbital 
Coulomb interactions, respectively, and Hj denotes the 
inter-orbital exchange interaction. Each term is repre- 
sented as 



We start with the Hamiltonian 



Hu ?^J7T«J7i ' 
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(10) 



and 



with 



Hj = jT.4aAb^'dja.'d,,. , (11) 
jaa' 

ijja = dj^^dj^^ and rij^ = J^'^'jicr- Here, we as- 
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sume that the three energy parameters are related as 
U = U' + J. By using the spin operator for the Cg elec- 
trons and the iso-spin operator describing the orbital de- 
grees of freedom, defined as 



(12) 
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Hon-^'/fe can be rewritten b" 



,yi 



Hn 



Coefficients of the spin and iso-spin operators, i.e., a and 
/?, are given by 



The bilinear terms with respect to the spin and iso-spin 
operators in the Hamiltonian is decoupled by introducing 
two kinds of auxiliary fields through the Stratonovich- 
Hubbard transformation. Then the partition function is 
rewritten as 

S= f Y[DSt,^M^,Dd,^,Dips^ipne- 1 ^^^"-'^''^^ (19) 
with 

Ld^^ da-yi (t) {dr - n) d„^i (r) 

£771 

+ XI ^77 ^^<T7^(T)rfa7'i(^") 
0-77' (ij) 

- 2aY,SeA^) ■ ^S^{t) - 2pY,f,{T) ■ 0T^{T) , (20) 



and 



(21) 



and 
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Being based on the above formulae (Eq.(|l^-(^l|)), we 
introduce the mean field approximation at this stage. At 
first, we consider the part of the Hamiltonian which de- 
scribes the t2g spin system, that is, the first and second 
terms in the right hand side in Eq. (|2l]). The mean field 
solution in this system is given by 



The minus sign in Eq. ( p^ ) means that the Coulomb inter- 
actions induce both spin and orbital (iso-spin) moments. 
In the path integral representation, the expression of the 
grand partition function is represented as 



S = y" Y[DSt,^,Dd,^,Dd,^, expj- J drLir) 



(16) 



with 



L{t)=H (t) + daJ^ (r) [dr - fl) d,^, (t) , (17) 



where r is the imaginary time introduced in the path in- 
tegral formalism, and dj^^-and d^^^ are the Grassmann 

variables corresponding to the operators dl^^ and d^^^, 
respectively. By utilizing Eq. (^3|) , Hamiltonian is rewrit- 
ten by 



(22) 



where (f>i is the solution of the following mean field equa- 
tion 



(23) 



By replacing the spin operator for the t2g spins by 
St2g,^y in Eq. (|l|) is given by 



(24) 
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Next, we focus on Ld in Eq. (|2C|). By using the momen- 
tum representation: 



da-f'j {t) 
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(33) 
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for X — S and T, we have 



(27) 



where is the Matsubara-frequency for fermion and the 
Gkk';nn';'ry';ai3 IS the Green function of the eg electron 
defined by 

'^fcfc';rm';77';a/3 = hi^n - fJ-) Snn' SapS^^' Skk' 
+ s'T' ^kk'Snn'SaP 



in terms of the doping concentration x. f (x) is the Fermi 
distribution function. 

By using the above expression of the free energy, we 
numericaUy calculate the spin and orbital phase diagram 
at zero temperature. We consider four kinds of spin align- 
ment in the cubic cell: F-, A-, C- and G-type. The pos- 
sibility of the spin canting is also discussed later. As for 
the orbital degrees, their ordering is represented by the 
alignment of the iso-spin. We specify the orbital state by 
the angle in the orbital space as follows: 



cos - \d^2_y2) + sin - |d3^2_ 



(34) 



^ ■ ifxik- fc', Un - u;„')Sai3 , (28) 



with 



1 \^ ,jV -ikRi+ik' Rj _ c'''^' X, , 



(29) 



After integrating over the Grassman variable, the parti- 
tion function is rewritten as 



i5{^}exp ^Tr In G^^.„„.^^.„^ 



drL^ 



(30) 



Then we adopt the mean field approximation by replac- 
ing the two kinds of auxiliary field, i.e., (ps and 0t by 
their values at the saddle point Lps and (pT- Finally, we 
obtain the expression for the free energy in the mean field 
approximation as 



which describes the direction of the iso-spin moment 

f = (- sin 0, 0, cos 61) . (35) 

We also consider four types of orbital ordering, i.e., F-, 
A-, C-, G-type, in the cubic cell. The angle in the or- 
bital space 6 is varied for each sublattice, and these are 
denoted as Oj and On in the / and // sublattices, respec- 
tively. Henceforth, we often use the notation such as, 
orbital G : (Sx^ - r"^) / {iy"^ - r^) = (G : 7r/3, -7r/3), 
through the relations. 



W3) I , 
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) + |c(3^2_r2) 

(37) 



Therefore, we consider the 4 (spin) x 4 (orbital) types of 
ordering with (0/, 0//), and numerically compare the free 
energy between them. 



i^ln [l + exp{-/3 (^^''^^-m)} 



(31) 



where E'^'^'> is the z^-th eigenvalue of Mkk'-ji'-ap defined 
by 



Mfcfc';77';a/3 = £fc^ ^kk'Saff Jt^^O'P ' Vsik - k')S^j> 

V N 



^=^17' ■ 0T{k- k')Sap ■ 



(32) 



Chemical potential fj, is determined by the following con- 
dition 



III. NUMERICAL RESULTS 
A. Parameters in the model 

The values of the energy parameters a, f3, Jh , Js^to, 
used in the numerical calculation are chosen as follows. 
In LaMnOa, to is estimated-by the photoemission exper- 
iments to be to 0.72eVj2J which we choose the unit 
of the energy below {to = 1). By employing U = 6.3eV 
and J — Jff = ]-OeV as those relevant to the actual 
manganese oxides parameters for the electron-electron 
interactions in the present model are a = 8.1, (3 = 6.67 
{a//3 — 1.21). The numerical calculation are also per- 
formed by using different sets of the enp|gw-.parameters 
in order to compare the previous works.E3~E3 There, the 
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effective Hamiltonians are derived by excluding the dou- 
bly occupied state in the Cg orbials. The superexchange 
interaction between nearest neighboring spin and orbital 
in these models are characterized by the energy in the in- 
termediate states in the perturbatioa. process. There are 
three kinds of intermediate stateSjCXESc^ i.e., the two 
electrons occupy 1) the different orbitals with the par- 
allel spin (the energy is U' — J), 2) the different ones 
with the antiparallel spins ([/' -f- J), and 3) the same or- 
bital ([/). The connection between these energies and 
the present energy parameters are roughly estimated as 
U' - J ^ (3, 5, and t/' + J - a + /3 from Eq.Q. 
Koshibae et alr3 discussed the orbital ordering at a; = 
by using the exact diagonalization method in the limit of 
U' — ^« U, U' + J corresponding to a/P » 1. Shiina 
et al.cn also studies the spin and orbital structure in the 
wide range of the parameters. In order to compare the 
above results, we show the two cases, that is, (case (A)) 
with a = 70 and /3 = 2.5 (a//3 » 1), and (case (B)) 
with a ~ 8.1 and (3 = 6.67. 



B. Undoped {x = 0) case 



Let us consider the undoped case. 



FIG. 2. 



In Fig.||, we show the calculated free energyF ( J5, x — 0) 
at X = as a function of J5 for each spin alignment in 
the case of (5//3 » 1 (case (A)). In each spin ahgn- 
ment, the orbital structure are optimized. In the case 
of Js = 0, the F-type spin alignment are favored due to 
the ferromagnetic superexchajage interaction under the 
doubly-degenerate orbital.cXEEl With increasing Js, the 
stable spin structure is changed from F- to A- , C- and 
finally G-type. This sequential change of the spin struc- 
ture is consistent with tlae|ptiieoretical studies based on 
the effective HamiltoniantllEa in the hmit a/ (3 » 1. 

FIG. 3. 

Spin and orbital structures are depend on Js and these 
are depicted in Fig.^. Rearrangement of the orbital strufe 
ture with changing Js, which is previously pointed out,EHI 
are also found. 

It is worth noting that the orbital structures also de- 
pend on the ratio a/p. 

FIG. 4. 

In Fig.^, the angle in iso-spin space {0i,9n) obtained 
in the spin A phase is presented. For a/ (3 » 1 the 
configuration {61, 9jj) = (90, —90) is obtained as the sta- 
ble orbital structure. In this orbital ordering, the en- 
ergy gain due to the super exchange process in the fer- 
romagnetic bonds, i.e., I {U' — J) takes its maximum. 



Here, t represents the effective electron transfer inten- 
sity in this superexchange process including the effects of 
orbital. On the other hand for d//3 « 1, the superex- 
change processes in the AF bonds characterized by t'^/U 
and jiU' -I- J), become important. Actually (Qi, 9u) — 
(180,-180) is the most preferable configuration in this 



parameter region. Orbital G : l3x 



^)/(3y^ 



((6*/, 6'//) = (60, — 60)), which is supposed experimentally, 
can not be the most stable for any a/f3. Rather than this 
structure, G:{y^ - z^)/{x^ - z^) ((6*/, 6i//)=(120, -120)) 
becomes stable around a/ (3 ^ 1.1. 

FIG. 5. 

These numerical results are understood by comparing 
the energy gains due to the superexchange processes be- 
tween the two orbital configurations shown in Fig. |^. 
In the processes with the energy of t^/{U' -f J) and 

/{U' — J), the transfer integral t between the occupied 
and the unoccupied orbital is concerned and it takes the 
same value for both configurations. On the other hand, 
in the process with t^/{U), relevant transfer is the one 
between the occupied orbitals along c-axis, which is al- 
ways larger for orbital: G:[y^ — z^) /(x^ — z^) than that 
for orbitahG : (3x^ — r^)/(3y^ — r^). Then there is no 
chance for orbital G: (3x^ — r^)/(3y^ — r^) to be the 
most stable structure for any a/ (3, when only the super 
exchange mechanism is considered. Hence we conclude 
that the JT coupling plays an indispensable role for the 
G:{ix^ — r^) / (3y^ — r^) orbital ordering at x = 0, which 
we will discuss in the next section. 

In the vicinity of a/ (3 ^ 1 , A- type spin ordering ap- 
pears even at Js = 0. This is consistent with the result 
obtained in the effective Hamiltonian at x = O.a In this 
region, the ferromagnetic interaction characterized by the 
energy /{U' — J) is suppressed, on the other hand, the 
AF ones with /U and / {U' + J) are enhanced. When 
a/ (3 approaches to unity with fixing Ju as several eV, the 
spin F structure at Jg = remains. It is concluded that 
the strong Jh plays an important role for the emergence 
of spin F at J5 = 0. 



C. Effects of the lattice distortion 



With considering the experimental fact that t 



the. static 

JT distortion rapidly disappears around x ~ 0.1,I1jEj we 
examine the effect of the JT distortion on the spin and or- 
bital phase diagram at x = 0. The JT distortion directly 
affects the orbital configuration. As a result, the phase 
boundary between spin F and spin AF phases is modified. 
To examine how the phase boundary is changed due to 
the JT distortion, we chose the parameters as a/ (3 » 1, 
i.e., case (A), where the phase boundary exists at x = 
as shown in Fig. 9. The JT coupling is expressed as the 
coupling between the iso-spin operator Tj and the local 
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lattice distortion Qj as follows:N 
HjT^9j2 -^^ + 111 ^oQ' + E ^(^j) ' (38) 

j 3 ] 

(39) 

where V{Qj) is the anharmonic potential for Qj. Instead 
of minimizing the total energy, we assume Qj as the ex- 
perimentally observed one. In a MnOe octahedron, Q^ 
is expressed as 



= Tj ( sin 0,(32.2„„2 + cos 0j (53^2 



(40) 



where Q^2_y2(^'^^2_^2^ is the base of the normal coordi- 
nate of the cubic- symmetric system defined as Qx^-y'^ = 
^ (Aj; - Ay) and Q3z2_^2 = {2Az - - Ay), and 
Aa denotes the elongation toward the a-direction (a = 
x,y,z). In this notation, the first term in Hjt, which is 
termed Hei-ph, is expressed as. 



el—ph 



\g\ 'S^ r,Vi ■ T, 



3 



3 ' 



with 



' sin Q 




(42) 



We choose the sign of the coupling constant -I- \g\ so 
that the d3^2_j.2 orbital may be stabilized for Q3^2_^2, 
consistent with the negative charge of oxygen ion. By 
the X-ray diffraction experiment, it is confirmed that 
the MnOe octahedron is elongated along the x- or y- 
direction and these octahedron are alternatively aligned 
in the a6-aplane.t3 In the present formula, it corresponds 
to r-j = r, = -60° for j e /-sublattice and 6^ = -60° 
for j e //-sublattice, i.e., (9/, 6//)= (60°, -60°). By 
adding H^i-ph to the Hamiltonian introduced in the pre- 
vious section (Eq.(|l|)), an additional molecular field for 
the iso-spin is introduced. As a result, the matrix ele- 
ment Mkh'-y-y'-afi in Eq.(|3^) is modified as 

— s'l'^ ^kk'SajS — Otac(3^ s^k~k' ^qs^ii' 



B Isl ^ Ba\ 



= -cos(-120°)sin6'/ 



,.BA 



.A 



: sin (-120°) cos 6*7 
: cos (-120°) cos 6*7 
= sin (-120°) sin 6I7 , 



(fiff — ifiT sm 



ip^ = ipT COS 



- tin 



II - C77 



II - an 



2 / W7 — W77 



Lpj^ = LfT COS 



and qs {qe) denotes the wave vector for the spin (or- 
bital) ordering. As well as the energy splitting between 
the two Cg orbital due to the JT distortion (Eq.(p])), 
the distortion modifies the transfer integral through the 
modification of the bond length I between Mn and O ions. 
According to the pseudo potential theory)^ the overlap 
integral between Mn 3c? and O 2p orbitals is proportional 
to f^^. Therefore, the variations of the transfer integral 
between Mn Sd orbitals is evaluated by using the param- 



(41) eter r(= rj) as 



i77 

i i+x(y) 



(r) = 



t'^"' (r 

ii-\-'x{y) V 



0) 



.^(l + 2r)^(l-r)^ ' 



and 
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(r) 



(1 



(44) 



(45) 



where we used the expression for the bond lengths in 
the elongated and shorten bonds as hong = lo {I + 2r) 
and Ishort = lo{l — r), respectively. We also consider the 
change in the magnitude of J5 due to the JT distortion. 
With the relation Jg (r) cx (r) , following relations are 
derived 



Js (r = 0) 



(l + 2r)^(l-r)' 



and 



Js ir) = 



Jsir = 0) 
(1 



\14 



(46) 



(47) 



where Jg (r) is the superexchange interaction along a di- 
rection with the distortion r. According to the X-ray 
diffraction experiment ,113 bond lengths are reported as 
hong = 2.14A and Ig^^ort — 1-98 or 1.96 A corresponding 
to r = 0.028 . In order to distinguish the two kinds of 
modification due to the JT distortion, that is, the en- 
ergy level splitting and the modification of the transfer 
intensity, we examine these effects separately by two pro- 
cedures as changing the value of g with fixing r and vice 
versa. Even for g = 0, the modification of the transfer 
(^"^Intensity lifts the degeneracy. 



where 



FIG. 
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In Fig.^ the stable orbital structure is shown as a func- 
tion of the diagonal coupling gr with fixing r. For (7 = 0, 
(6*/, Oil) is determined so as to lower the center of mass 
in the valence band. For sufficiently large g compared 
with to 7 JT distortion forces the orbital configuration to 
be (6'/,6'//)=(60, —60). Types of the stable orbital are al- 
most saturated as {ix^ — r^) /(3?/^ — r^) for gr/t^ > 0.5. 
When the Coulomb interactions are not introduced, this 
value gr/to ~ 0.5 is not enough to make the wave func- 
tions to be [3x^ — r^)/ (3?/^ — r^). Since the orbital is 
already polarized by the strong Coulomb interaction, the 
role of JT coupling is to rotate the direction of its po- 
larization. It is much easier than to induce the polar- 
ization. We concluded that the wave function is almost 
(3x^ — r^)/{3y^ — r^) aX,^^ 0, which in principle can be 
tested experimentally.EZlES 

FIG. 7. 

In Figj^, the variation of the phase boundary between 
spin F- and spin A-phases are presented as functions of 
r and gr. The value of the superexchange interaction at 
the boundary is termed Js{FA), hereafter. In the case of 
g = (Fig. ^ (a)), spin F is stabilized with increasing r. 
This is reasonable because the modification of the trans- 
fer intensity, described by Eqs. (^J) and (^5|), enhances 
the ferromagnetic superexchange interaction along the c- 
axis. On the other hand, introducing of the energy split- 
ting represented hygr with fixing r stabilizes the spin 
A structure (Fig. [t] (b)). In order to understand this 
results, we estimate the energy gain due to the energy 
splitting as follows. The energy difference between the 
spin F and spin A phases are represented as 

2S^JsiFA;g)^EA{g)-EF{g) (48) 

where Eji(p-){g) is the energy for the spin-F(A) with the 
JT coupling g, and Js {FA;g) is the superexchange in- 
teraction at the phase boundary. The prefactor 2 in the 
left hand side is came from difference of the number of 
the antiferromagnetic bonds between two phases. By us- 
ing this expression, the change of the phase boundary 
between the g = and g = cx) is estimated by 

2S^ (Js iFA;g - 0) - Jg {FA;g^ ex.)) 
= {EA{g = 0) - EAig = oo)) 
- {EF{g^O)~Epig^oo)) . (49) 

When the right-hand side of Eq.(|4^) is positive, the phase 
boundary Jg {FA;g) is increased with decreasing g. Be- 
cause the ferromagnetic superexchange interaction is only 
relevant at x = with the condition a//3 » 1, Ep(^A){9) 
is proportional to the sum of square of the transfer inten- 
sity between the nearest neighboring occupied and unoc- 
cupied orbitals (^0-11(5)), that is, 

i^F(A)(5)«- Cn(ff), (50) 

Ferro— bonds 



where ^ implies the summation over the ferro- 

Ferro— bonds 

magnetic bonds. 

FIG. 8. 

These quantities are tabulated in Fig||. As a result, we 
obtain 

{Js (FA; g = 0)-Js {FA; g = 00)) 
= 0.875 > , (51) 

with 

We conclude that Jg {FA; g) increases with decreasing 
g as shown in Fig.^ (b). This results implies that with 
relaxing the JT distortion, a frozen of the orbital config- 
uration is melted from (0/, Ou) = (60,-60), the spin F 
phase is stabilized in comparison with spin A phase. It is 
consistent with the experimental results where the spin- 
A phase is replaced by the spin-F insulator accompanied 
with reduction of the JT distortion as increasing x, al- 
though the present calculation is limited in the undoped 
case. 



D. Doped (2: 7^ 0) case 

In this subsection, we show the results in the finite hole 
doped case. 

FIG. 9. 



FIG. 10. 

In Fig.|| and |l^, we present the phase diagrams as a 
function of hole concentration x in the cases of (A) 
(a = 70,/3 = 2.5) and (B) (5 = 8.1, /3 = 6.67), respec- 
tively. Nonmonotonic behavior of the phase boundary is 
attributed to changes of the orbital structures. In both 
cases, the global features of phase diagram are almost 
the same. As discussed before, a/ (3 » 1 corresponds to 
U' — J « U, U' + J and the large Hund coupling. In this 
limit, the superexchange process for the AF interaction 
is neglected and the analyses of the calculated results be- 
come easier. Therefore, at first, we focus on the results 
in the case (A). 

It is found that in nearly the whole doping region ex- 
cept for X = 0, {3z^ — r^) and {x^ — y^) orbital struc- 
tures are stabilized in the spin A- and C-type phases. 
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respectively, because these orbitals are favorable to max- 
imize a gain of the kinetic energy in each spin struc- 
ture. Deviation from (x^ — j/^) structure in the A- 
type AF phase is found in 0.45 < a; < 0.75 and it is 
attributed to the hybridization between the occupied 
and unoccupied bands as discussed later. For spin G- 
type AF, the energy does not depend on the orbital 
so much, because the electron motion is blocked in 
all directions. In the spin F structure, on the other 
hand, the transfer is allowed in any direction. Or- 
bital structure in spin F changes continuously as x in- 
creases from orbital G: {x^ — y^^ /(3z^ — r^) near a; = 
to orbital F: {x^ — y^^ for x ^ 0.3, and to orbital 
A:{[3z^-r^] + [x^ ~ y^]) / {[Sz^ - r^] - [x^ - y^]) for 
0.3 < a; < 0.8 and finally orbital F: (Sz^ - r^). The or- 
bital structure is sensitively changed by changing x in 
comparison with that in the other spin structures. 

In order to understand the variation of the spin and 
orbital structures in the finite hole concentration region, 
let us consider the density of states (DOS) for each orbital 
configuration. 

FIG. 11. 

In Fig.|ll|, we present the schematic picture of DOS for 
several values of a and f3. In the case of large a and (3, 
the density of states is split into four bands. Each band 
are characterized by the direction of the spin and iso- 
spin, that is, the spin is parallel or anti-parallel and the 
iso-spin is parallel or anti-parallel to their mean fields, 
respectively. It accommodates an electron per site, and 
the lowest band corresponds to the state where both 
spin and iso-spin are parallel. Energy difference be- 
tween the two bands where the (iso)spin are parallel and 
anti-parallel, respectively, is given by a(ps= d{l — x)/2 
{Pi^T— P — x)/2). In the small doped case where 
(1 ~ x) » tt^l P is satisfied, the large energy gap ap- 
pears and the hybridizations between the bands are neg- 
ligible. The lowest band is therefore constructed almost 
only from the orbital described by Eq.(^. On the other 
hand, in sufhciently large doped case (1 — a;) ~ io//3j 
the energy gap shrinks and the orbital polarization is 
reduced. This change of the band gap with varying the 
hole concentration is able to be detected as the inter-band 
transition in the optical measurements. 

Next, we demonstrate how the orbital structure con- 
trols the dimensionality of DOS. 

FIG. 12. 

In Fig.|l^, we present DOS calculated in the several or- 
bital structures. The spin structure is assumed as the F- 
type and the hole concentration is fixed as a; = 0. At first, 
we focus on the case for orbital G:{x'^ — y"^) / (3z^ — r^) 
and F:(x'^ — y"^) (Fig|l| (b)) which are reahzed in x < 0.3 
in the ferromagnetic region as shown in Fig.0. Re- 



sults in the both cases are essentially the same. DOS 
shows a two-dimensional character, because of absence 
of the hopping integral along the c-axis. In this case 
there is no difference in the kinetic energy between spin 
F and spin A, and Jg favors spin A. This causes a 
dip structure in the phase boundary between spin F 
and spin A at x ~ 0.3 and leads an remarkable dif- 
ference from the prediction by the conventional double 
exchange model. In Fig.na (c), DOS for the orbital 
A : ( [3z2 - r2] + [x^ - y^]JJ{ [3z^ - r^] - [x^ - y^\ ) is 
shown. This orbital structure is stabilized in 0.3 < a; < 
0.8 as shown in Fig.^ Although DOS should be es- 
sentially three dimensional one if there is no hybridiza- 
tion, results has two peaks which resembles that in the 
quasi one-dimensional system (Fig.|2| (a)). This seems 
to be originated from the hybridization with the un- 
occupied band. For each case in Fig.[l^, the width in 
the lowest band is the same, as expected in the case of 
(1 — a;) » IqIP- Therefore, by adjusting the orbital 
structure, the shape of DOS is modified and the cen- 
ter of mass for the occupied states are changed so as to 
minimize the kinetic energy. From this viewpoint, a di- 
mensionality of the lowest energy band plays an essential 
role on a gain of the kinetic energy. In the regions of 
a: 1, one-dimensional-like dispersion is advantageous 
as shown in Fig.^. 

Let us consider the case (B) in FigJTo| where the more 
realistic energy parameters are adopted. At the moment, 
a value of Jg can not be estimated accurately, but there 
are two rough estimates. One is from tha,Neel tempera- 
ture Tat = 130X for CaMnOg (a: = 1.0), til which suggest 
Jg = T/v/7.5 1.7meV ^ 0.0023to in the mean field 
approximation. The fluctuations lower T/y, and hence in- 
crease the estimate for J5. Another estimate is obtained 
from the numerical calculations for LaMaOa (a; = 0.0), 
which suggests Js = 8meV = O.Ollio-Eil Although Js 
might depend on x in real materials, we tentatively fix 
Js to be 0.009to represented by the broken line in Figjic|. 
Then the spin structure is changed as A — > F — > A^ C 
^ G, as a; increases, which is in good agreement with the 
experiments shown in Fig.^c). 

As we mentioned above, the ferromagnetic phase is 
roughly divided into the two regions: the low doped 
region (a; < 0.3) and high hole doped one (a: > 0.3). 
The former is not reproduced by the calculation without 
the Coulomb interaction between Cg electrons, on other 
hand, the latter is not changed. We conclude that the 
superexchange interaction discussed in the previous sub- 
section and the conventional double exchange interaction 
are dominant in the lower and higher doped regions, re- 
spectively. In the present formulation, character of the 
superexchange interaction in the metallic phase is de- 
rived by the following mechanism. We consider the state 
where the AF spin structure or the AF orbital struc- 
ture exist. In this case, the eigen energy of Mkk' ■,'y-y' -,0113 
in Eq.(|3^) corresponding to the Hartree-Fock energy is 

roughly expressed as Ek ~ J e{kY + U^ff^xY , where 
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e{k) and Ueff{x) are the diagonal and ofF-diagonal ma- 
trix elements of Mkk>;-^-y'-ap , respectively. The former 
is the order of the transfer intensity and the latter is 
roughly estimated as Ueff{x) — Ucf ~ C^c(l — x). Lp is 
the auxiliary field for the spin and/or iso-spin degrees 
and Uc is the order of U . Therefore, the band width of 
the band is given by uj = E^^^^ — Ueff{x). ko is the mo- 
mentum where e{k) = 0. It is approximated as /Uc for 
small X. This results is in contrast to the large x case 
where the band width is the order of tj-. As de Gennes 
have pointed out in the case for small xp the kinetic en- 
ergy is determined by the product of the band width and 
the carrier concentration expressed in the present case as 
AE ^ wx — {t^ /Uc)x . We stress that the energy scale 

/Uc corresponds to that in the superexchange interac- 
tion. As increasing x, this ' superexchange character ' 
in the interaction is gradually replaced by the double ex- 
change one. Actually, the peak in the phase boundary at 
x ~ 0.15 in Fig.|l^ grows up with increasing t/Uc (Fig.p^ 
(a) ), on the other hand, the structure for x > 0.3 is 
almost unchanged. 

The modification of the F- and A-spin phases by 
changing the energy parameters discussed above exalains 
the recent experiment in (Lai-2Ndz)i_a;Sr2:Mn03.E3 The 
phase transition between the spin A and spin F metallic- 
phases were studied by changing the band width, which 
is controlled by z, and the hole concentration x. With 
increasing the band width, the critical hole concentration 
xf-a, where the phase transition occurs, is increased. 

FIG. 13. 

This experimental results is consistent with the present 
calculation shown in Figjl^ (b) where and Js oc t"^ are 
changed with fixing a and (3. It is found that the criti- 
cal carrier concentration xf-a is shifted to the higher x 
region with increasing to as consistent with the experi- 
ments. 

As shown in Fig.^ and the orbital structure in the 
spin F phase is sensitively changed with the energy pa- 
rameters and the carrier concentration. This implies that 
there are many nearly degenerate orbital configurations 
in this phase. 

FIG. 14. 

In order to investigate the situation in detail, we compare 
the free energies with assuming several orbital structures 
in the ferromagnetic phase (Fig.jlJ (a)). We vary the an- 
gle in the orbital space with assuming the ferromagnetic 
orbital configuration. The calculation is also performed 
in the A-type AF spin case (Figjl| (b)). It is found that 
the energy variation is an order of magnitude smaller in 
the case of spin F compared with the case of spin A. 
The difference between two cases are interpreted from 
the view point of the anisotropy of the electron trans- 
fer under the orbital ordering as follows. In spin A case, 



(x^— y^) orbital is realized to maximize the kinetic energy 
gain and the hopping along c-axis is forbidden. In spin 
F case, on the other hand, such kind of lowering of the 
dimensionality does not occur because the three crystal- 
lographic axes are equivalent in this spin structure. As 
a result, the orbital configurations have more freedom. 
This is the same physical idea behind th£_orbital liq- 
uid scenario proposed by present authors O where the 
two-dimensional orbital fluctuation, characterized by the 
(x^ — ?/^), (y^ — z^) and (z^ — x^) orbital ahgnments, 
is suggested. Through the orbital fluctuation among 
them, the kinetic energy in every direction is lowered 
and the ferromagnetic phase is stabilized. As a results, it 
is thought that the dip structure in the phase boundary 
between F- and A-spin structures shown in Fig.^ dis- 
appear and the two ferromagnetic phases in the low and 
high doped regions are connected. 

Another possibility, which enhances the ferromag^ 
netism in the region of a; < 0.175 is the JT distortion.Ea 
In experiments, however, the static JT distortion diSi 
appears rapidly with increasing of x around x ~ 0.1. Eil 
Then the dynamical JT distortion should be considered 
in the metallic-spin F state. According to the study in a 
large-c? modelJiJ the large Coulomb iiiteraction is essen- 
tial to explain both the_isotope effecica and the Raman 
scattering experimented Then the dominant role of the 
Coulomb interaction assumed in this paper seems to be 
reasonable. 

One of the most remarkable results in Fig.^ is an 
emergence of the spin A metallic phase for 0.2 < a; < 
0.5 and the spin C phase for a; > 0.5. These two 
phases are found in experiments using high quality sam- 
ples with narrowciiJband width: Pr1_2.Sr2.MnO3 and 
Ndi_2:Sra,Mn03.&By Also in Lai_2;Sr2;Mn03, spm A 
metallic is recently found to emerge for x > 0.54.E3 It 
is worth noting that the metallic spin-A phase is highly 
contrast to the spin-A insulating phase observed around 
X ^ 0. An existence of the spin canting in the metallic 
spin-A phase is theoretically excluded as follows. Accord- 
ing to the study by de GennesJJ the spin canting from the 
insulating A-type spin structure lowers the energy by 



TP —TP ^^inter 



Sz\Js\S^' 



(52) 



where Ecant is the energy in the spin canted phase, binter 
is the hopping integral along the c-axis, z is the coordi- 
nate number along this axis, and N is the number of the 
ion. The canting angle O is given by 





cos — 
2 



(53) 



From the consideration, the spin A phase around < a; < 
0.1 in Fig.|l^ is replaced by the spin canted phase and a 
value of Js at the phase boundary is corrected downward 
as 



Js [F-Cant) = Js {FA) 



36z\Js\S^ 



(54) 
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On the other hand, in the spin A metalhc phase (0.2 < 
X < 0.45) with the orbital structure (x^ — y^), the hop- 
ping along the c-axis is forbidden by the orbital struc- 
ture and hence no spin canting occurs. This is consis- 
tent with the recent neutron scattering experipipnt show- 
ing no canting in this spin A metallic phase.O Further- 
more, the large anisotropy in the jxesistivityt2l and the 
distortion of the MnOg octahedrono observed in the A- 
type AF at x=0.60 and the C-type AF at x=0.75 in 
Ndi_a:Sr2;Mn03 are consistent with the calculated or- 
bital structure: (a;^ — y^) for spin A and (3z^ — r^) in 
spin C. 

IV. DISCUSSION AND CONCLUSIONS 

Here we discuss the relation between the present work 
and the njjevious ones performed_a.t x = 0.0. Kugel and 
Khomskiio and Koshibae et aLE3 have dealt with the 
spin and the orbital orderings at x = 0, taking into ac- 
count the orbital degeneracy and the electron-electron 
interactions. They used the effective Hamiltonian ob- 
tained by the second-order perturbative expansion in the 
limit of strong Coulomb repulsion. Kugel and Khom- 
skii studied the ground state spin and orbital structures 
in the small limit of J/U' without the antiferromag- 
netic interaction J5 between t2g spins. In the mean 
field calculation, the A-type spin structure is reproduced 
but the orbital structure is almost (z^ — x'^)/{y'^ — z^) 
type ordering, which disagrees with one expected from a 
type of the observed JT distortion. On the other hand, 
Koshibae et al. also studied the spin and orbital struc- 
tures in the large limit of J/U' with taking into account 
Js- By using the exact diagonalization method in the 
finite cluster system, it is found that the spin correla- 
tion changes as F ^ A ^ C ^ G-type as Js in- 
creases. They also found that in spin A phase, component 
of (3:e2 _ ^2)^(3^2 _ ^2) (^2 _ x^Yl^y^ _ ^2^ Q^^ital 

alignments are enhanced, although both components are 
not distinguished in the orbital correlation function cal- 
culated there. Our calculation covers both of these two 
cases with the following features; a) Js is taken into ac- 
count, b) perturbative expansion is not used, i.e., ap- 
plicable for any parameters, a/ (3, J/U'. c) spin and 
the orbital orderings in the infinite system are studied. 
As shown in Figs. 1-3, in the limit of a/(3 » 1, a se- 
quential change of the spin structure with changing J5 
is observed. The orbital structure is also rearranged 
and only in the spin A-AF phase, {9i,9n) = (90°, 270°) 
is stabilized. On the hand, in the limit of a/ (3 ^ 1, 
{61,911) = (120°, -120°) is stabilized in the spin A-AF 
phase. Therefore, the results obtained in the previous 
calculationcjEj are reproduced by the present calcula- 
tion in the unified fashion. In addition, our calculation 
shows that the orbital ordering expected from a type of 
JT distortion can not be realized for any value of a/(3. 
It implies indispensability of the anharmonicity from the 



JT distortiorJla for proper description at a; = 0. Ttis is 
in accordance with the results-by the Hartree-FockEH and 
the first-principle calculations except that spin A phase 
is realized even without JT shown as in Fig.|lO|. 

Next, let us pay our attention to the doped case. As 
mentioned above, the perfectly polarization of the or- 
bital moment derived by the electron-electron interac- 
tions plays an essential role on the spin ordering. It con- 
trols the dinemsionality of the conduction band through 
the anisotropy of the transfer intensity. If not, the 
anisotropy of the conduction band is weakened, because 
of the large hybridyzation between the lowest band and 
the other ones. 

FIG. 15. 

As a comparison we show in Fig.^ a phase diagram 
without the orbital polarization by assuming (3 = 0. 
The phase diagram is dominated by the ferromagnetic 
state for reasonable values of Jg, and the nonmonotonic 
behaviors shown in Figs.^ and |l^ disappear. We con- 
clude that the almost saturated orbital polarization is 
essential to obtain the experimentally observed phase di- 
agram and the unique character observed in the each 
spin phase, e.g., the two dimensional conduction and the 
no spin canting in the metallic A-AF phase. For such a 
large orbital polarization in the metallic phase, the strong 
Coulomb interaction is indispensable rather than the JT 
coupling. 

In summary, we have studied the phase diagram of 
Ri_a:Aa;Mn03 at zero temperature in the plane of x (hole 
concentration) and Js(AF exchange interaction between 
t2g spins) in the mean field approximation. The global 
features are understood in terms of the superexchange in- 
teraction and the double exchange interaction, which is 
considerably modified from the conventional one due to 
the strong correlation and the orbital degeneracy. The 
large orbital polarization originated from the electron- 
electron interaction is indispensable to reproduce the 
phase diagram experimentally observed. 
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Figure captions 

Figure |]. Phase diagrams of (a) Lai_a;StJVIn03,i (b) 
Pri_^Sr^Mn03,El and (c) Ndi_:rSr:rMn03.El 

Figure ^. Free energies for each spin ahgnment as a 
function of the antiferromagnetic interaction between t2g 
spins J5 at a; = 0. The energy parameters are chosen to 
be a = 70 and f] — 2.5 (case (A)). Orbital structures are 
also noted. 

Figure ||. Orbital structures at x = as a function of 
the antiferromagnetic interaction between the spins 
Js- The energy parameters are chosen to be a = 70 and 
P = 2.5 (case (A)). 

Figure |4[ Orbital structures in each sublattice at a; = as 
a function of the electron-electron interaction parameters 
a/ p. The A- type spin alignment and the G-type orbital 
one are assumed. 

Figure p[ Orbital orderings described as 

(e»^,6'//)=(120°,-120°) and (60°, -60°) respectively. 

Figure ^. Orbital structures in the each sublattice as 
a function of magnitude of the JT splitting gr/to- The 
energy parameters are chosen to be a = 70 and /3 = 2.5 
(case (A)), r is fixed to be 0.028. 

Figure |^. Variation of the phase boundary between spin 
F- and spin A-phases as a function of (a) the magnitude 
of the JT distortion r with fixing the JT splitting to be 
g — (left-hand side panel), and (b) the JT splitting 
gr/to with fixing the magnitude of the distortion to be 
r — 0.028 (right-hand side panel). 



Figure |[ Magnitude of the transfer integrals for each 
orbital orderings, appearing in the calculation of Eq.(51). 



Figure |[ Mean field phase diagram as a function of the 
carrier concentration (x) and the antiferromagnetic in- 
teraction between <2g spins {Js )■ The energy parame- 
ters are chosen to be 5 = 70 and /3 = 2.5 (case (A)). 
Schematic orbital structure in the each phase is also 
shown. 



in each band, (a) a = /3 = case, (b) a 7^ 0, /? = case, 
and (c) a ^ 0, (3 ^ case. 

Figure Density of states of the lowest band at x = 
0. The orbital structures are assumed to be (a) 6 = 
180° (d3^2_^2) , (b) 6* = 0° {d^2_y2), and (c) = 90°. 
The ferromagnetic spin structure is assumed. The energy 
parameters are chosen to be a = 70 and (3 = 2.5 (case 
(A)). A Value beside each bond represents the transfer 
intensity. 



Figure [13[ Variation of the phase boundary Js{FA) be- 
tween the spin F and A-phases. (a) Values of the param- 
eters a and f3 are changed, (b) Values of the parameter to 
is changed. The cross point xfa represents the transition 
from Spin F to Spin A, where Js scales as Js oc t^ 

Figure |l^. The energy as a function of the orbital state 
characterized by 6 in the several value of x. (a) Spin F 
is assumed, (b) Spin A is assumed. In both cases, the 
orbital F-type structure is assumed. The energy param- 
eters are chosen to be a = 8.1 and /3 = 6.67 (case (B)). 

Figure |l^. Mean field phase diagram calculated with 
assuming /3 = 0. 



Figure |l^. Mean field phase diagram as a function of the 
carrier concentration (x) and the antiferromagnetic inter- 
action between t2g spins ( Js). The energy parameters are 
chosen to be a = 8.1 and /3 — 6.67 (case (B)). Schematic 
orbital structure in the each phase is also shown. Dotted 
line {Js =0.009) well reproduces the change of the spin 
structure experimentally observed (see text). 

Figure ^ Splitting of the band structure for eg elec- 
trons. The up and down arrows represent the direction 
of spin, while \a) and \b) represent one of the eg orbitals 
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